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ABSTRACT 


Independent Earth orientation measurements taken by the space-geodetic techniques of lu- 
nar and satellite laser ranging, very long baseline interferometry, and the Global Positioning 
System have been combined using a Kalman filter. The resulting combined Earth orientation 
series, SPACE2014, consists of values and uncertainties for Universal Time, polar motion, 
and their rates that span from September 28, 1976, to April 14, 2015, at daily intervals and is 
available in versions with epochs given at either midnight or noon. The space-geodetic mea- 
surements used to generate SPACE2014 have then been combined with optical astrometric 
measurements to form two additional combined Earth orientation series: (1) COMB2014, 
consisting of values and uncertainties for Universal Time, polar motion, and their rates that 
span from January 20, 1962, to April 14, 2015, at daily intervals and which are also available 
in versions with epochs given at either midnight or noon; and (2) POLE2014, consisting of 
values and uncertainties for polar motion and its rate that span from January 20, 1900, to 
March 23, 2015, at 30.4375-day intervals. 
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INTRODUCTION 


Reference series of Earth orientation parameters (EOPs) obtained by combining indepen- 
dent measurements of the Earth’s orientation are generated annually at the Jet Propulsion 
Laboratory (JPL) in support of tracking and navigation of interplanetary spacecraft. This 
report describes the generation of the most recent such combined Earth orientation series: 
SPACE2014, COMB2014, and POLE2014. Since the procedures used to generate these most 
recent series are similar to those used to generate previous such combinations, only a brief 
description of their generation is given here. Further details regarding the approach used at 
JPL to annually combine independent measurements of the Earth’s orientation can be found 
in Gross (1996, 2000) and Gross et al. (1998) 
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SPACE2014 


Data Sets Combined to Form SPACE2014 

SPACE2014 is a combination of independent space-geodetic measurements of the Earth’s ori- 
entation. Table 1 lists the space-geodetic series used in generating SPACE2014, giving their 
identifiers, the number of measurements from each series that were actually incorporated 
into SPACE2014, and the time interval spanned by those measurements. Note that the Uni- 
versity of Texas Center for Space Research (UTCSR) satellite laser ranging (SLR) Universal 
Time (LIT) values were not used in generating SPACE2014 due to problems associated with 
separating this component of the Earth’s orientation from the effects of unmodeled forces 
acting on the satellite that cause the node of its orbit to drift (see Gross et ah, 1998, p. 217) 
for further discussion about this point). For similar reasons, the International Laser Ranging 
Service (ILRS) satellite laser ranging lcngth-of-day (LOD) values have also not been used in 
generating SPACE2014. 

Since it was desirable to combine only independent measurements of the Earth’s orien- 
tation, only those series listed in Table 1 were used, even though other space-geodetic series 
are available from other analysis centers. When more than one series determined by the 
same measurement technique was used, care was taken to make sure that the measurements 
themselves were not included more than once. In particular, polar motion measurements 
from the JPL Global Positioning System (GPS) series were only used until the start of the 
International Global Navigation Satellite Systems (GNSS) Service (IGS) combined series 
EOP(IGS) 95 P 01 on January 1, 1995; polar motion measurements from the IGS com- 
bined series EOP(IGS) 95 P 01 were then used until the start of the IGS combined series 
EOP(IGS) 95 P 02 on June 30, 1996; polar motion measurements from the IGS combined 
series EOP(IGS) 95 P 02 were then used until the start of the accumulated IGS Solu- 
tion Independent Exchange (SINEX) combined series EOP(IGS) 00 P 03 on February 26, 
2000; and polar motion measurements from the accumulated IGS SINEX combined series 
EOP(IGS) 00 P 03 were used thereafter. Similarly, LOD measurements from the IGS com- 
bined series EOP(IGS) 95 P 02 were used until it ended on February 25, 2000 with LOD 
measurements from the accumulated IGS SINEX combined series EOP(IGS) 00 P 03 used 
thereafter. In addition, LIT measurements from the National Oceanic and Atmospheric Ad- 
ministration (NOAA) International Radio Interferometric Surveying (IRIS) Intensive UT1 
series were used until it ended on December 31, 1994; measurements from the Linked States 
Naval Observatory (USNO) National Earth Orientation Service (NEOS) Intensive LIT1 series 
were then used until it ended on December 4, 2000; and measurements from the Goddard 
Spaceflight Center (GSFC) NEOS Intensive LIT1 series were used thereafter. Finally, polar 
motion measurements from the LITCSR SLR series EOP(CSR) 96 L 01 were used until it 
ended on February 3, 1996, with measurements from the ILRS combined SLR series being 
used thereafter. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The Earth orientation series listed in Table 1 were first preprocessed by removing leap seconds 
from the Universal Time (UT1) values and, when necessary, by correcting the LIT1 values 
to be consistent with the extended definition of Greenwich Sidereal Time (GST) as adopted 


2 



Table 1 . Data Sets Combined to Form SPACE2014* 


Data Set Name 

Data Analysis 
Type Center 

Reference 

Data Span 

Number 

LLR (JPL15M01; VOL, UTO) 
McDonald Cluster 

LLR 

JPL 

Williams et al. (2012) 

Oct. 5. 1976, to Jan. 5, 2012 

747 

OCA 

LLR 

JPL 

Williams et al. (2012) 

Apr. 7, 1984, to Apr. 13, 2015 

1286 

Haleakala 

LLR 

JPL 

Williams et al. (2012) 

Feb. 10, 1985, to Aug. 11, 1990 

70 

Apache Point 

LLR 

JPL 

Williams et al. (2012) 

Jun. 4, 2006, to Nov. 17, 2014 

430 

UTCSR (CSR96L01) 
LAGEOS Polar Motion 

SLR 

UTCSR 

Eanes & Watkins (1996) 

Sep. 27, 1976, to Feb. 3, 1996 

2220 

ILRSA (08Mayl5; Polar motion) 
ILRS Primary Combination 

SLR ASI-CGS 

Sciaretta et al. (2010) 

Feb. 4, 1996, to Apr. 13, 2015 

6970 

DSN (JPL97R01; T, V) 
California-Spain Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Nov. 26, 1979, to Apr. 10, 2015 

1111 

California-Australia Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Oct. 28, 1978, to Sep. 30, 1997 

698 

NOAA (NOAA95R02) 
IRIS Intensive UT1 

VLBI 

NOAA 

Ray et al. (1995) 

Apr. 2, 1984, to Dec. 31, 1994 

2395 

USNO (N9903) 

NEOS Intensive UT1 

VLBI 

USNO 

Eubanks et al. (1999) 

Jan. 4, 1995, to Dec. 4, 2000 

1497 

NASA/GSFC (GSF2014a) 
NEOS Intensive UT1 

VLBI 

GSFC 

NASA/GSFC, VLBI Group (2015) 

Dec. 6, 2000, to Apr. 13, 2015 

4525 

NASA/GSFC (GSF2014a) 
Multibaseline 

VLBI 

GSFC 

NASA/GSFC, VLBI Group (2015) Apr. 12, 1980, to Apr. 14, 2015 

4895 

NASA/GSFC (GSFC1122) 
Westfort-Fort Davis 

VLBI 

GSFC 

Gordon et al. (1999) 

Jun. 25, 1981, to Jan. 1, 1984 

105 

Westfort-Mojave 

VLBI 

GSFC 

Gordon et al. (1999) 

Mar. 21, 1985, to Aug. 6. 1990 

18 

GPS (21APR04; Polar motion) 
Post-processed Flinn Analysis 

GPS 

JPL 

Heflin et al. (2004) 

Jun. 10, 1992, to Dec. 31, 1994 

817 

GPS (IGS95P01; Polar motion) 
IGS Final Combined (ITRF93) 

GPS 

NRCan 

Kouba & Mireault (1997) 

Jan. 1, 1995, to Jun. 29, 1996 

546 

GPS (IGS95P02; Polar motion) 
IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

Jun. 30, 1996, to Feb. 25, 2000 

1337 

GPS (IGS00P03; Polar motion) 
IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to Apr. 13, 2015 

5517 

GPS (IGS95P02; LOD) 

IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

Feb. 23, 1997, to Feb. 25, 2000 

1098 

GPS (IGS00P03; LOD) 

IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to Apr. 13, 2015 

5523 


* LLR, lunar laser ranging; JPL, Jet Propulsion Laboratory; VOL, variation of latitude; UT, Universal Time; OCA, Observatoire 
de la Cote d’Azur; UTCSR, University of Texas Center for Space Research; LAGEOS, Laser Geodynamic Satellites; SLR, satellite 
laser ranging; ILRS, International Laser Ranging Service; ASI, Agenzia Spaziale Italiana; CGS, Centro di Geodesia Spaziale; 
DSN, Deep Space Network; T, transverse; V, vertical; VLBI, very long baseline interferometry; NOAA, National Oceanic and 
Atmospheric Administration; IRIS, International Radio Interferometric Surveying; USNO, United States Naval Observatory; 
NEOS, National Earth Orientation Service; NASA, National Aeronautics and Space Administration; GSFC, Goddard Space 
Flight Center; GPS, Global Positioning System; IGS, International Global Navigation Satellite System (GNSS) Service; ITRF, 
International Terrestrial Reference Frame; NRCan, Natural Resources Canada; CODE, Center for Orbit Determination in 
Europe; IGN, Institut National de l’Information Geographique et Forestiere; SINEX, solution independent exchange; LOD, 
length of day. 
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by the International Earth Rotation and Reference Systems Service (IERS; IERS 1997, p. 
149). Since most of the series listed in Table 1 were already consistent with the extended 
definition of GST, this correction needed to be applied to only the NOAA IRIS Intensive 
UT1 series. 

Rotational variations caused by solid-Earth and ocean tides were also removed from the 
UT1 values. The effect of the solid-Earth tides was removed by using the model of Yoder 
et al. (1981); the model of Kantha et al. (1998) was used to remove the effect upon UT1 
of the ocean tides at the Mf, Mf', and Mm tidal frequencies. Since the Yoder et al. (1981) 
model already includes a contribution from the equilibrium ocean tides, only the Kantha 
et al. (1998) oceanic corrections to the Yoder et al. (1981) model were actually removed. 
Also note that the Kantha et al. (1998) model was used to remove the effect of ocean tides 
on only UT1, not on polar motion. Ocean-tide-induced polar motion variations were not 
removed from any of the polar motion observations. Finally, the only Earth orientation 
series listed in Table 1 that includes the effects of semidiurnal and diurnal ocean tides on 
the Earth’s orientation is the NOAA IRIS Intensive UT1 series. This series included these 
effects by adding to the released UT1 values the model of Herring (1993, also see Herring & 
Dong 1994). Hence, the same Herring (1993) model was used to remove them. 

On June 30, 1996, the IGS reference frame changed from the International Terrestrial 
Reference Frame ITRF93 to ITRF94; on March 1, 1998, it changed from ITRF94 to ITRF96; 
on August 1, 1999, it changed from ITRF96 to ITRF97; on February 27, 2000, it changed 
from ITRF97 to the IGS realization of ITRF97 known as IGS97; on December 2, 2001, it 
changed from IGS97 to the IGS realization of ITRF2000 known as IGS00; on January 11, 
2004, it changed from IGS00 to a new IGS realization of ITRF2000 known as IGbOO; on 
November 5, 2006 it changed from IGbOO to the IGS realization of ITRF2005 known as 
IGS05 (Ferland, 2006); on April 17, 2011 it changed from IGS05 to the IGS realization of 
ITRF2008 known as IGS08 (Rebischung et ah, 2011, 2012); and on October 7, 2012 IGS08 
was replaced by the new IGS realization of ITRF2008 known as IGb08 (Rebischung, 2012). 
These changes in reference frames potentially introduce discontinuities into the IGS combined 
Earth orientation series. However, the IGS Final combined series EOP(IGS) 95 P 02 and 
the accumulated IGS SINEX combined series EOP(IGS) 00 P 03 used here have had each 
of these segments aligned to the same IGS reference frame. Thus, to within the uncertainty 
in determining the corrections required to align each segment, the IGS combined series 
EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 used here should be reasonably consistent with 
each other. They were, therefore, concatenated with one common set of bias-rate corrections 
being determined for them, as described below. 

Adjustments Made to Space-Geodetic Series Prior to Combination 

Prior to combining the series listed in Table 1 to form SPACE2014, series-specific corrections 
were applied for bias and rate, the stated uncertainties of the measurements were adjusted 
by multiplying them by series-specific scale factors, and outlying data points were deleted. 
Values for the bias-rate corrections and uncertainty scale factors were determined by an 
iterative procedure wherein each series was compared, in turn, to a combination of all others. 
In order to minimize interpolation error (see Gross et ah, 1998, pp. 223-225), the comparison 
of each series to its reference combination was done at the epochs of the measurements of 
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that series by generating its reference combination using a Kalman filter that interpolates to 
(and prints its EOP estimates at) the exact epochs of those measurements. Also, both the 
bias-rate corrections and the uncertainty scale factors for all components of a given series 
were determined simultaneously in a multivariate approach using nonlinear weighted least 
squares. Using a multivariate approach allows the correlations between the components to be 
taken into account when determining the bias-rate corrections and uncertainty scale factors 
(see Gross et ah, 1998, pp. 22). 

All the series listed in Table 1 were included in the iterative procedure. Details of the 
iterative procedure are described in Gross (1996, 2000) and Gross et al. (1998) and will not 
be repeated here. These details include: 

1. Use of a reference series, SPACE2013 (Ratcliff & Gross, 2015), for initial bias-rate 
alignment; 

2. Analysis of each data type in its natural reference frame; 

3. Clustering of the McDonald hhR stations and, separately, the DSN very long baseline 
interferometry (VLBI) stations in California, Australia, and Spain; 

4. Initial convergence on values for the series-specific bias-rate corrections and uncertainty 
scale factors prior to the start of outlier detection and deletion; and 

5. Final convergence on these values after detecting and deleting all data outliers. 

At the end of the iterative procedure, relative bias-rate corrections have been determined 
that make the series agree with each other in bias and rate; uncertainty scale factors have 
been determined that make the residual of each series have a reduced chi-square near one 
when differenced with a combination of all others; and outlying data points (those with 
residual values greater than four times their adjusted uncertainties) have been deleted. A 
total of 94 data points, or about 0.22% of those combined, were thus deleted. 

For GSFC NEOS Intensive and GSFC Multibaseline VLBI, the initial alignment to 
SPACE2013 was performed individually for each baseline containing at least 90 data points 
and spanning more than 1 year. Baselines containing less than 1 year of data or fewer than 
90 data points were aligned as a group to SPACE2013. Subsequent steps of the iterative 
procedure were then conducted independently on the individual baseline groups. 

Starting April 17, 2011 the IGS changed their procedure for estimating uncertainty in the 
GPS EOP values. This introduced a change in the magnitude of the polar motion X (PMX), 
polar motion Y (PMY), and LOD uncertainties of approximately a factor of 6. To account 
for this sudden change, uncertainty scale factor values after April 17, 2011 were estimated 
separately from those of earlier epochs using a trial reference series. These trial scale factors 
were applied to the original IGS GPS data, which were then used in the iterative procedure 
described above to form a final reference series. The resulting scale factors shown in Table 2 
for the IGS00P03 data set include the contribution of both the trial scale factors and the 
overall iteratively determined scale factor. 
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Table 2. Adjustments to Space-Geodetic Series 3, 


Data Set Name 

Bias, 

mas 

LLR (JPL15M01) 
McDonald Cluster 
OCA 
Haleakala 
Apache Point 

VOL 

-0.054 ± 0.081 
-0.017 ± 0.032 
-1.804 ± 1.151 
-0.408 ± 0.250 

UT0 

-0.091 ± 0.086 
-0.014 ± 0.028 
-1.763 ± 0.692 
0.077 ± 0.264 

DSN (JPL97R01) 
California-Spain Cluster 
California-Australia Cluster 

T 

-0.810 ± 0.024 
0.582 ± 0.018 

V 

-0.139 ± 0.059 
0.364 ± 0.053 

NASA GSFC (1122) 
Westford-Fort Davis 
Westford-Mo j ave 

T 

8.858 ± 3.838 
0.629 ± 1.097 

V 

3.171 ± 5.918 
0.176 ± 1.675 


GSFC (GSF2014a) b 
G3GcKuRi 
HrRiWfWz 
MoRiWfWz 
Other Networks 

PMX 

0.074 ± 0.103 
0.775 ± 0.172 
0.203 ± 0.120 
0.022 ± 0.005 

PMY 

0.053 ± 0.103 
-0.206 ± 0.178 
-0.212 ± 0.115 
-0.069 ± 0.005 

UT1 

0.074 ± 0.168 
0.041 ± 0.284 
-0.065 ± 0.145 
0.077 ± 0.009 

NO A A (95R02) 
IRIS Intensive 

PMX 

PMY 

UT1 

0.070 ± 0.023 

USNO (N9903) 
NEOS Intensive 

PMX 

PMY 

UT1 

0.070 ± 0.023 

GSFC (GSF2014a) b 
KkNy 
KkWz 
TsWz 
NyTsWz 
KkSvWz 
Other Networks 

PMX 

PMY 

UT1 

0.721 ± 0.498 
0.025 ± 0.021 
-0.009 ± 0.057 
0.528 ± 0.180 
0.032 ± 0.121 
0.557 ± 0.209 

UTCSR (96L01) 
LAGEOS 

PMX 

-0.356 ± 0.010 

PMY 

0.740 ± 0.009 

UT1 

ILRSA (08Mayl5) 
Primary Comb. 

PMX 

0.042 ± 0.005 

PMY 

-0.070 ± 0.005 

LOD 

GPS (21APR04) 
JPL Post-Flinn 

PMX 

-0.201 ± 0.016 

PMY 

0.134 ± 0.013 

LOD 

GPS (IGS95P01) 
Final Combined 

PMX 

-0.026 ± 0.055 

PMY 

0.424 ± 0.054 

LOD 

GPS (IGS95P02) 
Final Combined 

PMX 

0.076 ± 0.004 

PMY 

-0.009 ± 0.004 

LOD 

-0.048 ± 0.010 

GPS (IGSOOP03) 
before 4/17/2011 
after 4/17/2011 

PMX 

0.076 ± 0.004 
0.076 ± 0.004 

PMY 

-0.009 ± 0.004 
-0.009 ± 0.004 

LOD 

-0.048 ± 0.010 
-0.048 ± 0.010 


Rate, mas/yr Uncertainty Scale Factor 


VOL 

UT0 

VOL 

UT0 

0.029 ± 0.019 0.007 ± 0.018 

1.398 

1.325 

0.000 ± 0.003 0.001 ± 0.003 

2.150 

1.789 

-0.373 ± 0.224 -0.392 ± 0.145 

2.777 

2.348 

0.024 ± 0.015 -0.004 ± 0.015 

1.728 

1.725 

T 


V 


T 

V 

0.187 ± 0.002 0.164 ±0.004 

2.123 

1.494 

0.086 ± 0.008 -0.064 ± 0.024 

1.437 

1.167 

T 


V 


T 

V 

0.859 ± 0.386 0.280 ± 0.592 

1.246 

0.910 

0.001 0.000 

2.501 

0.720 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

-0.152 ± 0.062 

0.120 ± 0.064 

-0.020 ± 0.102 

1.358 

1.369 

1.894 

0.074 ± 0.029 - 

-0.098 ± 0.030 

-0.072 ± 0.048 

1.175 

1.245 

1.704 

-0.067 ± 0.051 - 

-0.098 ± 0.048 

-0.071 ± 0.066 

1.432 

1.353 

0.460 

-0.001 ± 0.000 

0.001 ± 0.000 

-0.004 ± 0.001 

1.602 

1.580 

1.648 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 



-0.051 ± 0.007 



1.028 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 



-0.051 ± 0.007 



1.361 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 



-0.032 ± 0.026 



1.287 



-0.000 ± 0.001 



1.530 



0.006 ± 0.004 



1.776 



-0.026 ± 0.010 



1.805 



-0.001 ± 0.007 



1.736 



-0.029 ± 0.011 



1.580 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

0.095 ± 0.004 

0.122 ± 0.003 


0.773 

0.715 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.003 ± 0.000 - 

-0.000 ± 0.000 


3.144 

3.040 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.036 ±0.015 - 

-0.019 ± 0.012 


1.294 

0.951 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

0.169 ± 0.019 

0.245 ± 0.019 


1.258 

1.017 


PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.003 ± 0.000 

0.000 ± 0.000 

0.004 ± 0.001 

2.037 

1.654 

1.820 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

-0.003 ± 0.000 

0.000 ± 0.000 

0.004 ± 0.001 

2.037 

1.654 

1.820 

-0.003 ± 0.000 

0.000 ± 0.000 

0.004 ± 0.001 

11.713 

6.452 

11.111 


a) Reference date for bias-rate adjustment is 1993.0. See Table 1 footnotes, mas, milliarcseconds; PMX, polar motion X; PMY, polar motion Y. 

b) GSFC VLBI Station IDs from ftp://cddis.gsfc.nasa.gov/vlbi/ivscontrol/ns-codes.txt 

G3 = NRAO 85-3; Gc = Gilmore Creek; Hr = HRAS at Fort Davis; Kk = Kokee Park; Ku = 9m at Kokee Park; Mo = 12m at Goldstone; Ny = Ny 
Alesund, Svalbard, Norway; Ri = Richmond, FL; Sv — Svetloe, Russia; Ts = Tsukuba, Japan; Wf = Westford, MA; Wz = Wettzell, Germany. 
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Finally, each series was placed within a particular IERS reference frame by applying one 
final bias-rate correction that is common to all the series. The relative bias-rate corrections 
and uncertainty scale factors (determined by the iterative procedure described above) and the 
additional bias-rate corrections that align each series with an initial reference were applied to 
each individual series. These series were then combined and the resulting combination was 
compared to EOP(IERS) 08 C 04 (Bizouard & Gambis, 2011) during the interval 1984-2014 
in order to obtain the final corrections required to make it agree in bias and rate with the 
IERS series. 

The total bias-rate corrections and uncertainty scale factors determined by the procedures 
outlined above are given in Table 2. The values of the bias-rate corrections given in Table 2 
are the sum of: (1) the initial corrections applied to each series to align them to a common 
reference, (2) all of the incremental corrections applied during the iterative procedure, (3) 
the additional, common correction applied to the combined series in order to align the 
combined series with the IERS reference frame. The values of the uncertainty scale factors 
given in Table 2 are the products of all the incremental scale factors determined during the 
iterative procedure. The uncertainties of the bias-rate corrections given in Table 2 are the la 
standard errors of the incremental bias-rate corrections determined during the last iteration. 
There are no bias-rate entries in Table 2 for components that were either not used (e.g., the 
UTCSR SLR LIT1 component) or not available (e.g., the NOAA IRIS Intensive polar motion 
components). Note that the same common IERS rate correction is applied to all the data 
sets, including those such as the Westford-Mojave single baseline VLBI series, for which no 
relative rate correction could be determined. Therefore, the rate corrections given in Table 2 
for those series for which no relative rate corrections could be determined are simply the 
common IERS rate correction, but given in the natural reference frame of that series. In 
these cases, uncertainties for the rate corrections are not given. Also note that the entries 
for the bias-rate corrections in Table 2 for the IGS combined series EOP(IGS) 95 P 02 and 
EOP(IGS) 00 P 03 are the same. Since these entries were initially given within the same 
IGS reference frame, they were merged, and a common bias-rate correction was determined 
for them. 

Combined EOP Series: SPACE2014 

A Kalman filter was used to combine the series listed in Table 1 after the bias-rate corrections 
and uncertainty scale factors listed in Table 2 had been applied to them. The resulting 
combined Earth orientation series, SPACE2014, consists of values (Figure 1) and la standard 
errors (Figure 2) for polar motion, Universal Time, and their rates spanning September 28, 
1976, to April 14, 2015, at daily intervals; and it is available in versions for which the epochs 
are given at either midnight or noon. Leap seconds have been restored to UT1, and the 
effects of the solid Earth and ocean tides have been added back to the UT1 values using the 
same models for these effects that were originally used to remove them from the raw series: 
namely, the Yoder et al. (1981) model for the solid Earth tides and the Kantha et al. (1998) 
model for the ocean tides. However, semidiurnal and diurnal ocean tidal terms have not 
been added to and are therefore not included in the SPACE2014 UT1 values. 
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A Combined Earth Orientation Series: SPACE2014 



(smn) x uoijo^ .n^oj 



■xi" m fN ' — o 


(sui) Abq jo ipSuoq ssooxg 



o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

o 

o 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


CO 

<N 




<N 

co 

oo 

so 


(N 

<N 


so 

oo 


(sblli) x uoijoiai -re[ 0 <j (sui) mn- 1 in 


8 


Figure 1. Plots of the x-component of polar motion (la), y-component of polar motion (lb), UT1-UTC (lc), and excess 
length-of-day (Id) as given by the combined Earth orientation series SPACE2014. The discontinuous changes in the plot of 
UT1-UTC are caused by the presence of leap seconds. Note that the UT1-UTC values displayed in (lc) include the tidal 
variations, whereas the excess length-of-day values shown in (Id) do not. 




A Combined Earth Orientation Series: SPACE2014 
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Figure 2. Plots of the la formal errors of the x-component of polar motion (2a), y-component of polar motion (2b), UT1-UTC 
(2c), and excess length-of-day (2d) as given by the combined Earth orientation series SPACE2014. The inserts within panels 
(2a), (2b), and (2c) show that component’s post-1984 uncertainties on an expanded scale with the same units: milliarcseconds 
(mas) for polar motion, milliseconds (ms) for UT1-UTC. 






Figure 3 is a plot of the difference between the SPACE2014 polar motion, UT1, and 
LOD values and those of the IERS combined Earth orientation series EOP(IERS) 08 C 04. 
These two series are very consistent with each other, especially after January 1, 2000, when 
the root-mean-square (rms) of their difference is only 0.048 milliarcseconds (mas) for the 
x-component of polar motion, 0.043 mas for the y-component, 0.019 milliseconds (ms) for 
UT1, and 0.017 ms for LOD. Prior to 2000, the difference between these two series exhibits 
greater variability and even some systematic behavior, particularly in the x-component of 
polar motion. This systematic behavior is due to differences in the approaches used here and 
by the IERS to correct the bias and rate of the individual series before they are combined. 
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Difference of EOP(IERS)08C04 with SPACE2014 
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Figure 3. Plots of the difference between the IERS combined Earth orientation series EOP(IERS) 08 C 04 and SPACE2014 
formed by subtracting the SPACE2014 values from those of the IERS series. The difference between the x-component of polar 
motion is shown in (3a), the difference between the y-component is shown in (3b), the difference between UT1-UTC is shown 
in (3c), and the difference between the excess length-of-day is shown in (3d). 






COMB2014 


COMB2014 extends SPACE2014 by additionally incorporating the optical astrometric polar 
motion and UT1 series that was determined at the Bureau International de l’Heure (BIH) 
from an analysis of time and latitude observations by Li (1985, also see Li & Feissel 1986). 
This BIH optical astrometric series consists of values and uncertainties for polar motion and 
UT1 that span from January 5.0, 1962, to December 31.0, 1981, at 5-day intervals. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The BIH optical astrometric series was first preprocessed by removing leap seconds from the 
UT1 values and by correcting the UT1 values to be consistent with the extended definition 
of GST, as adopted by the IERS (IERS, 1997, p. 149). Rotational variations caused by 
solid Earth and ocean tides were also removed from the UT1 values. The same models that 
were used to remove the tidal effects from the series combined to form SPACE2014 were also 
used to remove them from the BIH series: namely, the Yoder et al. (1981) model for the 
solid Earth tides and the Kantha et al. (1998) model for the Mf, Mf', and Mm ocean tides. 
However, since the BIH UT1 measurements represent an average value over a 5-day-long 
observation window, and since 5 days is a substantial fraction of the monthly and shorter- 
period tides, the amplitudes of these solid Earth and ocean tidal terms were attenuated prior 
to their removal from the BIH UT1 measurements. (See Gross 1996, p. 8735 and Gross et al. 
1998, p. 226-227 for further discussion about this point.) 

Adjustments Made to BIH Series Prior to Combination 

The preprocessed BIH optical astrometric series was combined with the space-geodetic series 
that comprise SPACE2014 after Erst: 

1. correcting the BIH series to have the same bias, rate, annual terms, and semiannual 
terms as SPACE2014 

2. applying a constant multiplicative scale factor to the measurement uncertainties of 
the BIH series so that its residual, when differenced with SPACE2014, had a reduced 
chi-square of one; and 

3. deleting those data points, if any, for which residual values were greater than four times 
their adjusted uncertainties. 

The above adjustments were determined in a multivariate approach as was done for the 
series combined to form SPACE2014. The procedure used to determine these bias-rate 
and seasonal term corrections, uncertainty scale factors, and outlying data points has been 
described before (Gross, 1996, pp. 8735-8738) and will not be repeated here. The annual 
and semiannual terms of the BIH series were adjusted in order to correct for systematic, 
seasonally varying effects that are known to be present in optical astrometric measurements. 
Since the values of both the BIH series and the SPACE2014 series are given at midnight, 
interpolation error (see Gross et ah, 1998, pp. 223-225) is automatically minimized when 
differencing these two series for the purpose of determining the adjustments to be made to 
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the BIH series. Tables 3 and 4 give the resulting uncertainty scale factors and values and 
1 a standard errors of the corrections to the bias, rate, annual, and semiannual terms thus 
determined for the BIH series. When determining these uncertainty scale factors and the 
corrections to the bias, rate, and seasonal terms, no outlying data points were detected. 

Combined EOP Series: COMB2014 

A Kalman filter was used to combine the BIH series with the adjusted space-geodetic series 
that comprise SPACE2014 after first applying to the BIH series the uncertainty scale factors 
and corrections to the bias, rate, annual, and semiannual terms that are given in Tables 
3 and 4. The resulting combined Earth orientation series, COMB2014, consists of values 
(Figure 4) and la standard errors (Figure 5) for polar motion, Universal Time, and their 
rates that span from January 20, 1962, to April 14, 2015, at daily intervals and is available 
in versions with epochs given at either midnight or noon. Leap seconds have been restored 
to UT1, and the effects of the solid Earth and ocean tides have been added back to the UT1 
values using the same models for these effects that were originally used to remove them: 
namely, the Yoder et al. (1981) model for the solid Earth tides and the Kantha et al. (1998) 
model for the ocean tides. The full amplitude (i.e. , no tidal terms attenuated) of the effects 
of the solid Earth and ocean tides at the epoch of the time tag were added back to the UT1 
values. Semidiurnal and diurnal ocean tidal terms have not been added to and are, therefore, 
not included in the COMB2014 UT1 values. 


Table 3. Adjustments to Bias, Rate, and Stated Uncertainty of Optical Astrometric Series* 


Data Set 


Bias, mas 



Rate, mas/yr 


Uncertainty Scale Factor 


PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

BIH 

-2.426 ± 0.852 

-0.590 ± 0.699 

-7.825 ± 1.076 

1.051 ± 0.526 

1.182 ± 0.427 

5.592 ± 0.673 

1.810 

1.676 

1.943 

ILS 

-50.486 ± 2.357 

-2.956 ± 1.667 


0.427 ± 0.482 

-0.033 ± 0.340 


2.159 

1.502 



*Reference date for bias-rate adjustment of BIH series is 1980.0. Reference date for bias-rate adjustment of ILS series is 1970.0. 


Table 4. Adjustment to Annual and Semiannual Terms of Optical Astrometric Series* 


Data Set Coefficent of Sine Term, mas Coefficient of Cosine Term, mas 


PMX PMY UT1 PMX PMY UT1 


BIH 

Annual -5.528 ± 1.050 -7.215 ± 

Semiannual 2.282 ± 1.072 -0.506 ± 

ILS 

Annual -2.331 ± 3.288 8.447 ± 

Semiannual 0.240 ± 3.296 8.818 ± 


0.881 5.033 ± 1.387 -2.841 ± 

0.896 -0.532 ± 1.414 0.815 ± 

2.318 - — 8.103 ± 

2.325 - — 2.108 ± 


1.121 9.758 ± 0.922 -1.010 ± 1.443 

1.095 1.124 ± 0.903 1.941 ± 1.414 

3.324 -10.240 ± 2.347 

3.312 1.834 ± 2.337 


^Reference date for adjustment of annual and semiannual terms of BIH series is 1980.0. 
Reference date for adjustment of annual and semiannual terms of ILS series is 1970.0. 
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A Combined Earth Orientation Series: COMB2014 
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A Combined Earth Orientation Series: COMB2014 
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Figure 5. Plots of the la formal errors of the x-component of polar motion ( 5 a), y-component of polar motion ( 5 b), UT 1 -UTC 
( 5 c), and excess length-of-day ( 5 d) as given by the combined Earth orientation series COMB 2014 . 










POLE2014 


No optical astrometric observations taken at the stations of the International Latitude Ser- 
vice (ILS) were used when creating the BIH optical astrometric series used in COMB2014 
(Li, 1985; Li & Feissel, 1986). The ILS polar motion measurements (Yumi & Yokoyama, 
1980), which are based solely upon latitude observations made at the ILS stations are, 
therefore, independent of those comprising COMB2014 and have therefore been combined 
with them to form POLE2014. Being based solely upon latitude observations, the ILS series 
contains no UT1 measurements, but consists solely of polar motion measurements that span 
1899.8-1979.0 at monthly intervals. Although no uncertainties are given with the individ- 
ual polar motion values, the precision with which they have been determined is estimated 
to be 10-20 mas (Yumi & Yokoyama, 1980, p. 27). An initial uncertainty of 15 mas was, 
therefore, assigned to each of the ILS polar motion values. Since this assigned measurement 
uncertainty will be adjusted later, its initial value is arbitrary as long as it is not zero and 
serves merely as an a priori estimate to be used in the series adjustment procedure described 
below. 

The ILS series was combined with COMB2014 to form POLE2014 after: 

1. correcting the ILS series to have the same bias, rate, annual, and semiannual terms as 
COMB2014; 

2. applying a constant multiplicative scale factor to the measurement uncertainties of 
the ILS series so that its residual, when differenced with COMB2014, had a reduced 
chi-square of one; and 

3. deleting those data points, if any, whose residual values were greater than four times 
their adjusted uncertainties. 

These adjustments were also determined using a multivariate approach for fitting a bias, 
rate, and these seasonal terms to the difference of the ILS series with COMB2014 during 
1962.0 to 1979.0. The measurement uncertainties of the ILS polar motion values were ad- 
justed by determining and applying a scale factor that made the residual of this fit have a 
reduced chi-square of one. During this procedure to determine uncertainty scale factors and 
bias, rate, and seasonal term corrections, no outlying ILS data points were deleted since no 
data points had residual values greater than four times their adjusted uncertainties. Tables 
3 and 4 (in the COMB2014 section) also give the resulting uncertainty scale factors and 
values and la standard errors of the corrections to the bias, rate, annual, and semiannual 
terms thus determined for the ILS series. 

A Kalman filter was then used to combine the ILS series with the adjusted BIH and space- 
geodetic series that comprise COMB2014, after applying to the ILS series the uncertainty 
scale factors and corrections to the bias, rate, annual, and semiannual terms that are given 
in Tables 3 and 4. The resulting combined Earth orientation series, POLE2014, consists of 
values (Figure 6a and 6b) and ler standard errors (Figure 6c and 6d) for polar motion and 
its rate that span from January 20, 1900, to March 23, 2015, at 30.4375-day intervals. 
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A Combined Earth Orientation Series: POLE2014 
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Figure 6. Plots of the x-component of polar motion (6a), the y-component of polar motion (6b), the la formal errors of the 
x-component of polar motion (6c), and the lcr formal errors of the y-component of polar motion (6d) as given by the combined 
polar motion series POLE2014. 






DISCUSSION 


The Kalman filter that was used here to combine Earth-orientation measurements contains a 
stochastic model of the process that is used to propagate the state vector and its associated 
state covariance matrix forward in time to the epoch of the next measurement. For polar 
motion excitation, the stochastic model includes a random-walk term with equal noise forcing 
both the x- and y-components of excitation and originally having a white-noise power spectral 
density of 246.4 mas 2 /day (Morabito et al., 1988). This level of polar motion excitation 
process noise in the Kalman filter was increased to 739.2 mas 2 /day in order to better match 
the observed spectrum of polar motion excitation. Increasing the excitation process noise 
reduces the level of smoothing applied to the polar motion components of the propagated 
state vector and increases the covariance associated with those components. Thus, the 
SPACE2014, COMB2014, and POLE2014 polar motion and polar-motion excitation values 
are not as heavily smoothed as were those of SPACE2003, COMB2003, POLE2003, and 
earlier combinations produced at JPL; and the uncertainties assigned to the SPACE2014, 
COMB2014, and POLE2014 polar motion and polar-motion excitation values are somewhat 
larger. 

The LOD components of the IGS Final combined EOP(IGS) 95 P 02 and EOP (IGS) 
00 P 03 series were incorporated into JPL’s combinations starting with SPACE2007 and 
COMB2007. In general, GPS LOD estimates are contaminated by orbital artifacts that 
must be removed prior to combination with other LOD or UT1 measurements (Chin et al., 
2009). When determining the IGS Final combined LOD series, however, the IGS applies 
constraints that have the effect of removing these orbital artifacts. Therefore the IGS Final 
combined LOD series was incorporated into the Kalman filter as a true LOD measurement 
whose uncertainties are assumed to be “white” (i.e., random, uncorrelated in time, and with 
a flat power spectral density). Including the daily IGS Final combined LOD series in JPL’s 
combinations helps to compensate for less frequent UT1 measurements and was found to 
improve the agreement of JPL’s combined LOD estimates with independent atmospheric 
and oceanic angular momentum values, especially at the highest frequencies. 

Since a Kalman filter has been used to generate SPACE2014, COMB2014, and POLE2014, 
the resulting polar motion and UT1 values are smoothed to a degree that depends upon both 
the spacing between the measurements being combined (which determines how far the state 
vector and state covariance matrix must be propagated forward in time) and the uncertain- 
ties that have been assigned to the measurements. Since improvements to the observing 
systems, both in the hardware and software and in the number of systems, have led to in- 
creasingly precise determinations of the Earth’s orientation, and since the time resolution 
of the measurements has generally increased in concert with the measurement precision, the 
degree of smoothing applied to the SPACE2014, COMB2014, and POLE2014 values is a 
function of time, with the earlier values being more heavily smoothed than the more recent 
values. 

Daily EOP values are reported in SPACE2014 since the NOAA IRIS and GSFC NEOS 
Intensive UT1 values are given at daily intervals, as are the GPS and ILRSA combined 
SLR values (although gaps exist in each of these data sets). However, prior to the start of 
these data sets, the measurements combined to form SPACE2014 are given less frequently; 
therefore, the Kalman filter used to combine these measurements also interpolates them in 
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order to produce a series of equally spaced values. In order to be consistent with SPACE2014, 
daily EOP values are also reported in COMB2014 even though the BIH optical astrometric 
series used in COMB2014 is given at 5-day intervals. Thus, SPACE2014, COMB2014, and 
POLE2014 are equally spaced series of smoothed, interpolated Earth-orientation parameters. 

The combined Earth-orientation series SPACE2014, COMB2014, and POLE2014 are 
available from JPL’s Geodynamics and Space Geodesy Group via anonymous ftp at: 
<ftp://keof.jpl.nasa.gov/keof/combinations/2014> and upon request from the au- 
thors: Todd.Ratcliff@jpl.nasa.gov or Richard.Gross@jpl.nasa.gov. 
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Acronyms and Terms 


ASI 

BIH 

CGS 

CODE 

COMB 


DSN 

EOP 

G3 

Gc 

GNSS 

GPS 

GSFC 

GST 

Hr 

HRAS 

ID 

IERS 

IGN 

IGS 

ILRS 

ILRSA 

ILS 

IRIS 

ITRF 

IVS 

JPL 

Kk 

Ku 

LAGEOS 

LLR 

LOD 


Agenzia Spaziale Italiana 

Bureau International de l’Heure 

Centro di Geodesia Spaziale 

Center for Orbit Determination in Europe 

extends the SPACE series by additionally incorporating BIH optical 
astrometric measurements of polar motion and UT1 

Deep Space Network 

Earth orientation parameter 

IVS station code for 85-3 antenna at NR AO Green Bank, WV 

IVS station code for Gilmore Creek, Fairbanks, AK 

Global Navigation Satellite System 

Global Positioning System 

Goddard Space Flight Center 

Greenwich Sidereal Time 

IVS station code for HRAS 

Harvard Radio Astronomy Station at Ft. Davis, TX 
Identification 

International Earth Rotation and Reference Systems Service 

Institut National dc l’lnformation Geographique et Forestiere 

International Global Navigation Satellite Systems (GNSS) Service 

International Laser Ranging Service 

ILRS Primary Combination Solution 

International Latitude Service 

International Radio Interferometric Surveying 

International Terrestrial Reference Frame 

International VLBI Service for Geodesy & Astrometry 

Jet Propulsion Laboratory 

IVS station code for Kokee Park, Kauai, HI 

IVS station code for 9-m antenna at Kokee Park, Kauai, HI 

Laser Geodynamics Satellites 
lunar laser ranging 
length-of-day 
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mas 

MO 

ms 

NASA 

NEOS 

NOAA 

NRAO 

NRCan 

Ny 

OCA 

PMX 

PMY 

POLE 


Ri 

rms 

SINEX 

SLR 

SPACE 

Sv 

T 

Ts 

USNO 

UT 

UTO 

UT1 


UTC 

UTCSR 

V 

VLBI 

VOL 

Wf 

Wz 


milliarcsecond 

I VS station code for 12-m antenna at Goldstone, CA 
millisecond 

National Aeronautics and Space Administration 
National Earth Orientation Service 
National Oceanic and Atmospheric Administration 
National Radio Astronomy Observatory 
Natural Resources Canada 

I VS station code for Ny Alesund, Svalbard, Norway 

Observatoire de la Cote d’Azur 

polar motion X 
polar motion Y 

extends the COMB series by additionally incorporating 1LS optical 
astrometric measurements of polar motion 

1VS station code for Richmond, FL 
root mean square 

Solution Independent Exchange 
satellite laser ranging 

a combination of independent space-geodetic measurements of the 

Earth’s orientation 

IVS station code for Svetloe, Russia 

transverse component of Earth orientation from single-baseline VLBI 
IVS station code for 32-m antenna at Tsukuba, Japan 

United States Naval Observatory 
Universal Time 

Universal Time as determined at a particular observatory 
The principle form of Universal Time; mean solar time at 0° longi- 
tude obtained by correcting LITO for observatory location and polar 
motion. 

Coordinated Universal Time 

LIniversity of Texas Center for Space Research 

vertical component of Earth orientation from single-baseline VLBI 
very long baseline interferometry 
variation of latitude 

IVS station code for Westford, MA 
IVS station code for Wettzell, Germany 
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